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ABSTRACT 

GJ 436b is a Neptune-size planet with 23.2 Earth masses in an elliptical 
orbit of period 2.64 days and eccentricity 0.16. With a typical tidal dissipation 
factor (Q' ~ 10 6 ) as that of a giant planet with convective envelope, its orbital 
circularization timescale under internal tidal dissipation is around 1 Gyr, at least 
two times less than the stellar age (> 3 Gyr). A plausible mechanism is that 
the eccentricity of GJ 436b is modulated by a planetary companion due to their 
mutual perturbation. Here we investigate this possibility from the dynamical 
viewpoint. A general method is given to predict the possible locations of the 
dynamically coupled companions, including in nearby/distance non-resonant or 
mean motion resonance orbits with the first planet. Applying the method to 
GJ 436 system, we find it is very unlikely that the eccentricity of GJ 436b is 
maintained at the present location by a nearby/distance companion through 
secular perturbation or mean motion resonance. In fact, in all these simulated 
cases, GJ 436b will undergo eccentricity damp and orbital decay, leaving the 
present location within the stellar age. However, these results do not rule out 
the possible existence of planet companions in nearby/distance orbits, although 
they are not able to maintain the eccentricity of GJ 436b. 

Subject headings: celestial mechanics, methods: analytical, methods: N-body 
simulations, stars: individual (GJ 436), planetary systems: dynamics 



1. Introduction 



The discovery of the first extrasolar planet around a pulsar' 1 ', which was quickly followed 
by the detection of first Jupiter-like planet around the star 51 Peg^, opened a new era for 
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planetary science. The planet around 51 Peg is known as a hot Jupiter: planets with masses 
comparable to Jupiter's, orbits typically within 0.1 AU and surface temperatures ~ 1000 
K®. To date, more than 310 planets have been discovered in about 250 planetary system^. 
With the improvement of detection precision and the use of various techniques, the minimum 
mass of discovered planets around main sequence stars is down to 3.3 Earth masses (M©), 
which is MOA-2007-BLG-192-L b around a star with mass of O.O6M [ 4 l 

Among the planets discovered so far, GJ 436b is the only transiting Neptune-mass 
planet orbiting an M-type star. It was first detected by radial velocity techniques®, with its 
orbital elements refined later by Maness et al.®. The transiting signals of GJ 436b were first 
discovered by Gillon et al.^ and followed by lots of work' 8-13 ]. These observations reveal GJ 
436b with a mass of 23.17M e and a radius of 4.22i? (see Table 1) . 



Table 1: The parameters of GJ 436 and its planet companion GJ 436 



Parameter 


Value 


Ref. 


M* (M ) 


u.±oz_ 012 


1 


R* {Rq) 


n 4fi4+o.oo9 


1 


stellar age (Gyr) 


6 +4 


1 


M p (M e ) 


23.17 ±0.079 


1 


R p (i? e ) 


A 99+0.09 
^• zz -0.10 


1 


M p sini(Mj) 


0.0713 ±0.006 


2 


a(AU) 


0.0285 


2 


P(days) 


2.64385 ± 0.00009 


2 


e 


0.16 ±0.019 


2 



Refs.: 1. iTorred (120071 ). 2. Maness et all (120071 ). 



One of the most interesting characteristics about GJ 436b is its significant eccentricity 
(0.16) in an orbit (2.64 days) very close to the host star. Assuming a tidal dissipation factor 
(Q' ~ 10 6 ) as that of a gas giant planet, its orbit circularization timescale under internal 
tidal dissipation is around 1 Gyr. On the other hand, the fiducial age of the host star is 
6l 5 Gyr' 11 !, and according to observation, GJ 436 has low rotation velocity and does not 
exhibit particular strong chromospheric activity nor photometric variability®, indicating 
an age > 3 Gyr. As the orbit is not circulated by planetary tide, either GJ 436b has 
Q' > 6 x 10 6 , or there is a planetary companion which induces a periodic modulation of its 
orbital eccentricity. 
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Considering that radial velocities of GJ 436 reveal a long-term trend, Maness et al. 
proposed the presence of a long-period (~ 25 yr) planet companion with mass ~ 0.27Mj 
(Jupiter mass) in an eccentric orbit (e ~ 0.2) ^ . Recently, Ribas et al. suggested that the 
observed radial velocities of the system are consistent with an additional small, super-Earth 
planet in the outer 2:1 mean-motion resonance with GJ 436b' 14 ]. Such a possibility was also 
studied from the dynamical viewpoint t 15 '. More recent inspection of transit data implies 
that GJ 436b is perturbed by another planet with mass < 12M e in a non-resonant orbit of 
12days (0.08 AU)[ 16 l 

In this paper, we exam the possibility of a nearby or distant undiscovered planet through 
dynamical considerations. The key issue here is to locate the undetected companion that 
can both excite and maintain the eccentricity of GJ 436b. First we present the dynamical 
restrictions on the possible locations of the companion in section 2. Then in section 3 we show 
the analytical and numerical results on the eccentricity modulation by a companion in nearby 
or distant non-resonance orbits, followed by the investigation of the possible companion in 
resonant orbits in section 4. Conclusions are presented in section 5. Although the method 
is derived in GJ 436 system, it can be applied to any other systems in similar situations. 

2. Models and restrictions for planet companion 

Consider a planetary system of a star with mass M* and two planets with masses Mi 
and Mi- For simplicity we assume the two planets are in coplanar orbits and employ a 
general coplanar three-body model. The orbital elements of the two planets are denoted by 
dj, ej, Aj, Wi, which are semi-major axis, eccentricity, mean longitude, longitude of periapsis 
of planet i(i = 1, 2), respectively. The index i is labelled so that Mi is the inner planet with 
ai < 0,2 throughout the evolution. For the present problem, as it is unlikely to have a planet 
inside GJ 436b to modulate its eccentricity, we suppose M 1 is GJ 436b and let M 2 denote GJ 
436c. Suppose GJ 436b is located in an initially circular orbit with the present semi-major 
axis, we will study the mass range and the location of the planet companion that excites 
the eccentricity of GJ 436b to the observed value (0.16) due to either secular or resonance 
interaction. In this section we present three conditions that restrict the extension of mass 
and orbital elements for the undiscovered companion from (1) Doppler technique, (2) Hill 
stability consideration and (3) tidal circulation timescale. 
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2.1. Observational restriction from Doppler technique 

Radial velocity technique detects the stellar wobbles of orbital motion in the presence 
of a planet M p . In terms of the orbital elements of M p , the stellar radial velocity is given 

as! 17 !: 

[2ttG\ 1/3 M p sim 1 

V P J (M* + M P )V3 V ^3^I ^ ' 

where P is the period of the planet orbit. As M p <C M*, the above equation can be simplified 

as 



lOdays/ \ 10M e J \M Q J y/\ 



e 



Due to the perturbation of stellar photosphere, the limit precision that the Doppler technique 
can achieve is around 3 m/s' 17 '. Assume the companion we are to locate is comparable or 
below this limit, an undetect planet with stellar radial velocity K < 3 m/s has maximum 
mass of, 

/ p \ 1/3 

M P sini<5.9M e ^-j (1 - e^ 3 - (3) 



Using the elements and physical parameters in Table 1, and assuming sinz = 0.9784 
by the comparison of GJ 436b's mass from radial velocity and transit technique, we get the 
possible location of companion in the period-mass (P2-M2) space in Fig.l. 



2.2. Hill stability condition 



Dynamical stability of the planet system is the necessary condition for the presence of a 
second planet. One of the practical stability is Hill stability, which requires the ordering of 
two planets unchanged during the history of evolution. A coplanar, non- resonant, two-planet 
system is Hill stable if the following condition is satisfied^ 18 ' : 

2M T ^ , , o4/3 M X M 2 



H = 



G 2 M 3 



c 2 h > 1 + 3 4 / 3 - 



M* 2/3 (M! + M2) 4 / 3 



+ ... 



(4) 



where .\/ 7 M- ■ \I, ■ M* is the total mass of the system, M = M 1 M 2 + M 2 M* + M*M U 
C and h are the total angular momentum and energy of the three-body system, respectively. 
Since M\,Mi M*, by omitting up to the second order terms in mass ratios of planets 
to the star in the expression of total angular momentum and energy, the left-hand side of 
equation (j3J) can be approximated as, 



H 



1 



(m 1 +M 2 y 



M 1 Ja 1 (l-e 2 1 ) + M 2 



a 2 (l -e 2 2 ) 



M l M 2 
ai a 2 



+ 



(5) 
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P 2 (days) p 2 ( da V s ) 



Fig. 1. — Permissible regions (shadow areas) of the planet companion in period-mass (P-i — M2) 
plane with eccentricity e<i = in panel (a) or e2 = 0.2 in panel (b) under three restrictions: (1) 
radial velocity K = 3 m/s or 5 m/s, assuming sini = 0.9784(blue and green dot-dash lines, below 
which the radial velocity is smaller); (2) Hill stability (black solid line denoted by 'HS', the right 
part is Hill stable); (3) Tidal circulation timescale r c i rc = 6Gyr, assuming Q' 2 = 10 6 (red dashed 
lines, the right part with r c j rc > 6 Gyr). 

Denote m = Mj/M*, 73 = y/a^l — efj, the above criterion for the two planets in Hill stable 
regions reads' 19 ' , 

0*1 + /^r V71 + m 2 ) 2 ( - + -) > 1 + 3 4 /Vi/^i + ^\ (6) 

The Hill stability region by Eq. (jHJ) is also plotted in Fig.l. 



2.3. Restriction from tidal circulation timescale 

A close-in planet produces tidal bulges on the stellar surface, causing energy dissipation 
on the star and angular momentum exchanges between the stellar spin and planetary orbital 
motion. Meanwhile the star also generates tidal dissipation on the planet, resulting in an 
eccentricity damping and orbital decay' 20,21 ] . In the ideal case that both the stellar and 
planetary spins are aligned with the orbit, the secular evolution rate of the eccentricity can 
be expressed as follows I 22 ~ 24 l ( 

e = g P + g* (7) 
_ /3 (e) + H /l(e )(^)], (8, 



9p,* 



81 ne 



2 Q> \M P * 



M*. 



R 
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where g p and g* are contributions from the planet and star, respectively, n is the mean 
velocity of orbital motion, M* tP , Q'^ p , R* )P , fi*,p are the masses, the effective tidal dissipation 
factors, radii and spin rates of the star and the planet, respectively. Some functions of 
eccentricity used here and later are: 

/ 3 (e) = (1 + f e 2 + f e 4 + ^e 6 )/(l - e 2 ) 13 / 2 , 
/ 4 (e) = (1 + fe 2 + |e 4 )/(l - e 2 ) 5 , 

/ 6 (e) = (1 + f e 2 + ff e 4 + f e 6 + f| e 8 + I ffee 10 )/(1 + 3e 2 + §e 4 ), W 
fr(e) = (1 + ffe 2 + 8e 4 + || e 6 + ff e 8 + T f| f e 10 )/(1 + 3e 2 + fe 4 ). 

For close-in planets with tidal dissipation factor Q' p < 10 6 , dissipation in planets dominates^ 25 ' . 
Omitting contribution from the star in Eq. (j7]) and assuming the planet spin has reached the 
synchronization equilibrium (Q p ~ n), the timescale of orbital circularization (r circ = e/e) 
induced by planetary tidal dissipation is given as^ 25 !, 



Me 

And the associate timescale of orbital decay (rjecay = a/a) in elliptical orbits is 

(l-e 2 )Me) _ 

rdccay 2e 2 / 6 (e) circ ' 1 j 

For GJ 436b, with the elements and physical parameters in Table 1 and assuming Q' = 10 6 , 
the circularization timescale of GJ 436b is 1.0 Gyr, around five times less than the fiducial 
stellar age(6 Gyr), and r decay = 19 Gyr for GJ 436b at the present location with e = 0.16. 
So some mechanisms are needed to maintain its eccentricity during the evolution. 



3. Planet Companions in Non- Resonant Orbits 

Secular perturbations between two planets in non-resonance orbits exchange their an- 
gular momentum, thus modulate their eccentricities, leaving their semi-major axes almost 
unchanged. When tidal dissipation is present on either of the planet, eccentricity modulation 
is effective only when the timescale of secular perturbation (r sec ) is significantly shorter than 
that of the circularization (r circ ). With an octopole Legendre expansion in ratio a = a^ja-i, 
Mardling derived the period of secular perturbation at the limit of ei 1, including the 
effect of general relativity^ 26 ' , 



1 



-Pict 



-3 



M2 



-i 



e 3 
-2 



i r -u Ml 



+ i4 



(12) 
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where Pi is the mass of first planet (Mi), e 2 = \/l — e|, 7 = 4ct~ 3 (72iai/c) 2 (M i! /M2) is the 
ratio of general relativity to companion perturbation on periapsis precession of Mi, with 
ri\ the mean motion of M\ and c the speed of light. For a companion in a nearby orbit, 
Tscc <C Tcirc is easily fulfilled, as 7 1 and a is moderate. For a giant companion in a distance 
orbit, 7 ^> 1, thus Eq. (fl2l) can be simplified as r sec Pi/3/(niai/c) 2 = 1.5 x 10 4 yrs, 
which is independent of the planet companion and is much shorter than the circularization 
timescales of both planets. 

Under secular perturbation of M 2 , the maximum eccentricity of Mi that can be achieved 
from an initial circular orbit (eio = 0) \^ 2& \ 

(13) 

This equation can be used to locate a global approximate region of the planet companion in 
either nearby or distance orbits, while general three-body simulations should be performed 
to give a precise location. 

Eq. (1131) provides a mechanism that can excite e\ to a moderate value at M\ > M 2 . In 
fact, from the expression of ei max , there exists a singularity at a cr i t when the denominate of 
Eq. ffTB"]) is zero so that ei max tends to infinity when Mi > M 2 . In reality, as ei max becomes 
large enough, the above approximation ei <C 1 is no more valid and we should resort to 
numerical simulations. Fig.2a shows the maximum eccentricity (ei max ) as a function of 
companion's period and mass (P 2 -M 2 ) derived from Eq. fflHj) with e 2 o = 0.2. To verify this, 
we perform some three-body simulations with Mi in initial circular orbits, and obtain ei max 
under the modulation of M 2 at different locations. Fig. 2b shows the results by three-body 
simulations. The critical locations (a cr i t ) when the singularity occurs in Eq. (1131) are also 
plotted in the curves with the asterisks, showing a roughly good agreement between the 
analytical and numerical results. According to Fig.2, it is possible to excite ei max = 0.16 in 
close orbits with M 2 > 10M e and e 2 = 0.2. For gas giant companion, only nearby orbit is 
possible. We investigate these situations in detail as follows. 

Nearby Orbits. We perform extensive three-body simulations, including general relativity 
effect, on the initial P 2 — e 2 plane, with a nearby companion mass of 5M e or 10M e . The 
results are shown in Fig.3a,b. The shadow regions in Fig.3a,b are the most possible locations 
of the companion that can maintain e\ = 0.16 by secular perturbation, combined with the 
three restrictions present in section 2. As two examples, we present the evolution of two 
orbits from the permissible regions under tidal dissipation (with illustrative Q[ = Q 2 = 100) 
in Fig. 3c, d. The eccentricity of GJ 436b can be excited and maintained (with a periodic 
modulation) to 0.16 only within 10 5 years. Considering the linear dependence of tidal force 
on Qt, (i=l,2), these simulations indicate that ei max ~ 0.16 can be maintained only for 1 



eimax = -^ae 2 e 2 ' 
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Fig. 2. — Maximum eccentricity of Mi that can be excited at the present location by a planet 
companion with initial eccentricity e 2 = 0.2. Panel (a): 3-D plot of ei max in the plane of P2 — M2 
from theoretical Eq. f)13[) . Panel (b): ei max obtained from numerical simulations of a general 
three-body model. The asterisks denote the singularities from Eq. f)13[) . 

Gyr provide Q\ = 10 6 . 

Distant Orbits. According to equation (1T31) and numerical simulation (Fig. 2b), ei max is small 
unless M2 is in a highly eccentric orbit. Fig. 4a plots the region in P2 — €2 plane that a 
companion can generate a maximum eccentricity ei max = 0.16. They are calculated from 
equation (|T3l) and confirmed by full 3-body simulations. Thus it is almost impossible for a 
companion in orbit of P2 ~ 1 yr to produce ei max = 0.16. The two possible locations of M 2 
suggested by Maness et al.' 6 ' are also investigated and plotted in Fig.4d, which shows they 
can only excite negligible eccentricities of Mi. 



4. Planet Companion in Resonant Orbits 

For a conservative coplanar two-planet system, the motion of two planets can be de- 
scribed by a Hamiltonian system with four degrees of freedom, which is non-integrable. 
However, near a generic (p + q)/p (q 7^ 0) mean- motion resonance, the degrees of freedom 
of the system are reduced from four to two by averaging technique' 27 ' . Below we will show 
that, the minimum initial eccentricity of unseen planet M 2 can be deduced approximately 
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Fig. 3. — Permissible regions and orbital evolution for the companion M2 in nearby orbits that 
can excite e\ = 0.16 at the present location. Panels (a) and (b) are permissible regions (shadow 
areas) in P2 — ei plane obtained from general three-body simulations for companion mass 5M^ and 
lOM®, respectively. The boundary lines of the three restrictions (K= 3 m/s or 5 m/s; Hill stability 
denoted by HS, and r c i rc = 6Gyr with Q' 2 = 10 6 ) are also plotted. Evolution of Pi, e\ for two orbits 
(denoted by filled circles) in panels (a) and (b) are shown in panels (c) and (d), respectively, with 
initial parameters shown in the labels. Tidal dissipation factors of two planets are taken Q' 2 = 100 
in panels (c) and (d), with density of M2 taken as P2 = 3 g cm -3 . 
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0.28 
0.24 
0.20 
0.16 
0.12 
0.08 
0.04 
0.0' 



p 2 (yrs) 



-d: 0.3Mj, 1yr, e 20 =0.85 
c2: 0.3MJ, 1yr, e 20 =0.9 



(C) 




0.36 
0.32 
0.28 
0.24 
0.20 
0.16 
0.12 
0.08 
0.04 
0.0! 



4.5x1 0' 
4.0x10"' 
3.5x10"' 
3.0x10"' 
2.5x10"' 



■b1: 0.1 Mj, 0.2yr, e 20 =0.55 
b2: 0.1 Mj, 0.2yr, e 20 =0.7 



(b) 




2.0x10 3 4.0x10 3 6.0x10 3 8.0x10 3 1.0x10 4 
Time (yrs) 



-d1:0.12Mj, 14yr, e 20 =0.6 
d2: 0.27M y 25yr, e 20 =0.2 



(d) 



5.0x1 3 1.0x10 4 1.5x10 4 2.0x1 4 
Time (yrs) 




1x10 



2x10 3x10 
Time (yrs) 



4x1 4 



Fig. 4. — Permissible region and orbital evolution for distance giant companion Mi that can excite 
ei = 0.16 at the present location. Panel (a): Minimum eccentricities of companion with different 
mass M2 = O.lMj, 0.3Mj, lMj from Eq. (|13p and are confirmed by general three-body simulations. 
Panels (b-d): Evolution of e\ of six orbits notated in panel (a) by three-body simulations, with 
mass and initial elements shown in the labels. Orbits in panel (d) are proposed by Maness et alJ 6 '. 
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from the conservation of total orbital energy and angular momentum, with the help of the 
averaged Hamiltonian function. 

Adopt the following planar canonical variables' 27 ' 28 ^ , 



Ai, J 1 = L 1 + s(I 1 + I 3 ), 

A 2 , J 2 = L 2 -(l + s)(/i + / 2 ), 

ai = (1 + s)A 2 - sAi - zux, h = Li(l - a /1 - ef ), 

cr 2 = (1 + s)X 2 - sAi - rz? 2 , I 2 = L 2 (l - a/1 - ef). 



(14) 



where Aj, tUj are the longitude of mean motion, longitude of periapsis of Mj (i = 1, 2) 
respectively, and 

s=p/q, Li = Mly/tJ~a~i, (i = 1,2) 

/ii = G(M* + Mj), Aff- u ' u ' 



(15) 



with G the gravitational constant. The Hamiltonian H of the system can be expressed as: 

H = H + H U (16) 
where the first term corresponds to the two-body contribution given by: 

(17) 

1=1 4 

The second term, is the disturbing function. Up to the first order in the masses, it has 
the following expression t 28 l : 

.U;.U> .U;.U, . , . 

^i = ~ G ^^ h M (xix 2 + y x y 2 + ziz 2 ) (18) 

where A = (r 2 + r 2 — 2rir 2 cos S) 1 ^ 2 , and for the coplanar three-body system, S = f\ — 
f 2 + vj\ — w 2 with fi the true anomaly of the orbits m« (i = 1,2). In terms of the elements 
on (114p . all periodic terms in the Hamiltonian ( jTT51) contain only three independent angular 
variables <7i,<t 2 ,Ai — A 2 , thus the system is three degrees of freedom^ 27 !. The canonical 
moment conjugate Ai + A 2 is an integral of motion, namely J\ + J 2 = const. By averaging 
the synodic angle Q = (Ai — A 2 )/g, we obtain an averaged system with Hamiltonian function, 

H = H + ^- [ E x dQ (19) 



277^0 

In practice, the above averaged Hamiltonian can be obtained only numerically. 

The averaged system with Hamiltonian f|T9l) is of two degrees of freedom, with the energy 
being the only integral. To show all possible solutions in (ei, e 2 ) space for all possible phase 
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angles <7i, o 2 is impossible. So we fix only the symmetric resonance period orbits with initial 
Ci,^ being set at either or tc. Fig. 5a shows the energy level curves of Hamiltonian at 
10M e from equation ffTUj) on the (ei, e 2 ) plane. Based on the contour lines (or the averaged 
Hamiltonian function, equivalently) and fixed on symmetric period orbits, we derive the 
minimum initial eccentricity of the companion, e 2m i n , with which Mi can evolve to e± = 0.16, 
as a function of companion mass (Fig. 5b). Interestingly, e2 m in has a power-law dependence 
on the mass ratio Mi/M 2 , with an approximation relation: 

e 2min «0.14(^|) 1 / 2 (20) 

The relation holds for the three major resonances 2:1,3:1 and 5:2, and is independent of the 
location of M\. 




2E-7 6E-7 9E-7 1E-6 2E-6 2E-6 2E-6 3E-6 3E-6 



Fig. 5. — Panel (a): Constant-energy contour of the average Hamiltonian (fTUj) with M 2 = lOM®. 
Panel (b): Minimum eccentricity of M 2 that can excite e\ = 0.16 according to the energy level 
curves at 2:1,3:1 and 5:2 resonances. The red dashed curve is calculated by three-body simulations 
for the 2:1 resonance. 

The above relation can be deduced from the conservation of total angular momentum, 
J\ + -h = L\\J\ — ef + Li\f \ — e\ = const (independent of s). In fact, in an ideal situ- 
ation such that ei ma x (e2min) occurs at e 2 = 0(ei = 0, respectively), the conservation of 
total angular momentum requires, L 10 -\/l — ef max + L 2 q = L w + L 2 oa/1 — el min , which gives 
approximately, e 2min = y/L w / L 20 ei ms , x w (M 1 /M 2 ) 1/2 (n 2 /n 1 ) 1 / 6 e lmSLX , where L 10 ,L 20 are 
the corresponding elements evaluated at the resonance center. It shows that e 2min depends 



13 



weakly on the mean motion ratio, 
approximate formula fl20|) . 



At the 2:1 resonance, let e\ x 



0.16, we derive the 



From the above derivation of equation f[2"Ul) . we can see that, the specific resonance 
structure, which would be very complicated in a general three-body model, is not considered. 
So the relation (|20|) holds approximately only, and for a real e 2m j n , we shall resort to numerical 
simulations. Fig. 5b plots also the results from the three-body simulations including the 
general relativity effect. The discrepancy between the relation and the simulation is large 
especially when M 2 is small. 

Fig.6 shows the evolution of two typical orbits in 2:1 mean motion resonance with Mi 
for a companion M 2 = 5M e , as proposed by Ribas et alJ 14 '. The eccentricity they proposed 
is 0.2, below the value of e 2m in plotted in Fig.5b. As we can see, t\ can not be excited to 
0.16(Fig.6b). For the higher e 20 = 0.40 case, it can excite to e\ = 0.16 initially, but the 
eccentricity is damped soon. In both cases the orbit of GJ 436 decays to inner orbits. Other 
mean motion resonances show similar results, indicating that it is impossible to maintain 
ei = 0.16 by a resonant companion. 
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Fig. 6. — Evolution of Pi,e\ of two orbits initially in 2:1 resonance with Mi by three-body 
simulations. Companion mass M 2 = 5M^ and density p 2 = 3 g cm -3 , initial eccentricities of M 2 
are 0.2 and 0.4. The orbit with e 2 o = 0.2 was proposed by Ribas et alJ 14 '. Tidal dissipation factors 
of two planets are taken as Q' = 100. 
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5. Conclusions 

Perturbation from a companion M2 seems to be the most plausible approach to maintain 
the moderate eccentricity of GJ 436b in a close-in orbit. In this paper, we study extensively 
all different locations of the planet companion in the following three situations: 

(i) In nearby orbits. The eccentricity of GJ 436b can be excited to 0.16 with a broad 
range of planet mass (above few Earth- masses) . An interesting case is the second planet can 
have a mass much smaller than Mi. However, since the orbital decay time is short (~ 20Gyr) 
for GJ 436b at the present location with the observed eccentricity, significant orbital decay 
(~ 25%) is expect so that GJ 436b would be in a much closer orbit, and the eccentricity 
of GJ 436b would be damped within the stellar age. Thus the eccentricity of GJ 436b can 
not be maintained by a nearby companion. On the other hand, one can not rule out the 
possibility that companions exist in nearby orbits of GJ 436b as long as they are in stable 
orbits, e.g., the possible existence of a planet with < 12M e in an orbit of 12days from transit 
data' 16 l However, the companion, if exists, does not account for the eccentricity of GJ 436b. 

(ii) Distance Orbits. Distance companions on moderate eccentric orbits can not excite 
and maintain the significant eccentricity of GJ 436b. However, the presence of giant plan- 
ets in extended orbits (~ few AU) is also possible, although they can not modulate the 
eccentricity of close-in planets. For example, there may exist a long-period (~ 25 yr) planet 
companion with mass ~ 0.27 Mj inferred from radial velocities of GJ 436 ' 6 '. 

(Hi) Resonance Orbits. Although a border range of companion mass can excite the 
eccentricity of Mi to a moderate value, similar to situation (i), significant orbital decay 
would occur so that the two planets will leave 2:1 resonance soon, if GJ 436b has a normal 
dissipation factor {Q' ~ 10 6 ). In this case, the presence of a companion in the outer 2:1 
resonance with GJ 436b is unlikely. However, if we have to confess that GJ 436b has a 
extremely high dissipation factor (Q' > 6 x 10 6 ), the existence of companions in resonance 
orbits can not be ruled out in this extreme case, e.g., the planet (~ 5M e ) in the outer 2:1 
resonance place proposed by Ribas et aU 14 '. However, even it exists, it does not account for 
the eccentricity of GJ 436b. 

Based on the extensive investigations in this paper, we think ei = 0.16 of GJ 436b 
can not be maintained by a companion in either nearby or distance orbits through secular 
perturbation or mean motion resonance. Thus the maintaining of its eccentricity remains a 
challenge problem, unless GJ 436b has a extremely high dissipation factor {Q' > 6 x 10 6 ). 

Acknowledgments. We thanks Dr. Gregory Laughlin for useful discussions. 
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